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Abstract 

We study effects of noncommutativity on the phase space generated by a non-minimal scalar 
field which is conformally coupled to the background curvature in an isotropic and homogeneous 
FRW cosmology. These effects are considered in two cases, when the potential of scalar field has 
zero and nonzero constant values. The investigation is carried out by means of a comparative de- 
tailed analysis of mathematical features of the evolution of universe and the most probable universe 
wave functions in classically commutative and noncommutative frames and quantum counterparts. 
The influence of noncommutativity is explored by the two noncommutative parameters of space 
and momentum sectors with a relative focus on the role of the noncommutative parameter of 
momentum sector. The solutions are presented with some of their numerical diagrams, in the 
commutative and noncommutative scenarios, and their properties are compared. We find that 
impose of noncommutativity in the momentum sector causes more ability in tuning time solutions 
of variables in classical level, and has more probable states of universe in quantum level. We 
also demonstrate that special solutions in classical and allowed wave functions in quantum models 
impose bounds on the values of noncommutative parameters. 

PACS number: 04.20. - q ; 04.90. + e ; 04.20.Fy ; 98.80.Qc 
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1 Introduction 



Scalar fields are an integral part of modern models in particle physics [1], and recently play very 
important roles in cosmology and have become a powerful tool to build cosmological models as well. 
They have key role in some of these models as current models of early cosmological inflation [2], or, in 
the viability of scalar field models as favorite candidates for dark matter [3]. Scalar field cosmological 
models have extensively been studied in the literatures, see, e.g., Refs. [4] and references therein. In 
the simplest interactions, a scalar field is coupled to gravity. In many cosmological models, scalar 
fields present degrees of freedom and appear as dynamical variables of corresponding phase space, 
where this point can be regarded as relevance of noncommutativity in these models. 

The proposal of noncommutativity concept between space-time coordinates was introduced first 
by Snyder [5], and about twenty years ago, a mathematical theory, nowadays known as noncommu- 
tative geometry (NCG), has begun to take shape [6] based on this concept. In the last decade, 
study and investigation of physical theories in the noncommutative (NC) frame, like string and Mr- 
theory [7, 8], has caused a renewed interest on noncommutativity in the classical and quantum fields. 
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In particular, a novel interest has been developed in considering the NC classical and quantum cos- 
mology. In these studies, the influence of noncommutativity has been explored by the formulation of 
a version of NC cosmology in which a deformation of minisuperspace [9]-[12] or, of phase space [13] is 
required instead of space-time deformation. From qualitative point of view, noncommutativity in the 
configuration space leads to general effects, however, a non-trivial noncommutativity in momentum 
sector introduces distinct effects in what concern with the behavior of dynamical variables. 

Our purpose in this work is to build a NC scenario for the Friedmann-Robertson-Walker (FRW) 
cosmology including matter field via a deformation achieved by the Moyal product [7] in classical and 
quantum level. We introduce effects of noncommutativity by two parameters, namely and /?, which 
are the NC parameters corresponding to space and momentum sectors, respectively. Then, we will 
show that impose of noncommutativity in the momentum sector causes more ability in tuning time 
solutions of variables in classical level, and has more probable states of universe in quantum level. 

The work is organized as follows. In Section 2, we specify a model and inspect it in the classical 
version within the commutative and NC frames. Section 3 considers the quantum version of this 
model by investigating universe wave functions and compares their properties in the commutative 
and NC frames. A brief conclusion is presented in the last section. 



2 The Classical Model 

We consider a classical model consisted of a cosmological system that is presented by a four- 
dimensional action with a non-minimally coupled scalar field to gravity in a FRW universe. To 
specify the NC effects of the model, we first treat the commutative version, and then the NC one in 
the following subsections. 

2.1 Commutative Phase Space 

A general action for a non-minimally coupled scalar field can be described by 



A = f*=* 



d 4 x, (1) 



where g is the determinant of the metric g^, R is the Ricci scalar, V{4>) and /(</>) are potential and 
coupling functions of the scalar field, respectively. We assume a homogeneous scalar field, that is 
4> = 4>(t), and the following FRW metric of the minisuperspace 

ds 2 = ^2^2 + Q 2 {t) (^-^ + r 2 dfl^j , (2) 

where N(t) is a lapse function, a(t) is a scale factor and k specifies geometry of the universe. Substi- 
tuting the metric (2) in action (1), one obtains the Lagrangian density 

w(^4)-4 / +"Mj^). < 3 > 

where total time derivative terms have been neglected. 

We restrict our considerations to the case of conformally coupled scalar field model [12, 14]. The 
general reason for selecting such a scalar field is that it allows exact solutions in simple cases, as those 
discussed along this work, and it is rich enough to be useful as a probe for the significant modifications 
that NCG introduces in classical and quantum cosmologies. That is, we set /(</>) = 1/(2k) — £</> 2 /2, 
where k = SttG/c 4 and £ is the non-minimal coupling parameter that represents a direct coupling 
between the scalar field and curvature, and has an arbitrary value. Obviously, the case £ = is the 
minimally coupling situation, however, as mentioned, we consider the conformal coupling case, i.e. 
£ = 1/6, and employ the unites H = 1 = c and k = 3. 

Based on these assumptions and rescaling the scalar field as 
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X = a(f>/V2 



the Lagrangian (3) reads 



aa 2 av 2 kN\ 2 

C = kNa - — + -£ — - a NV . 

N N a 



Thus, the corresponding Hamiltonian is 



n = N (-^ + P A_ ka+ ^ + a Zy 

\ 4a 4a a 



(4) 



(5) 



where p a and p x are the canonical conjugate momenta. For the conformal time gauge selection, 
namely N = a, one gets 



n 



Then, the Hamilton equations are 



1 



(6) 



a = {a,H} = -^Pa, 
Pa = {Pa, H} = 2ka - 4a 3 V + V2 X a 2 V , 
X = {X,^} = -p x , 
P x = {p x ,H} = -2k X -V2a 3 V, 



(7) 



where the prime denotes derivative with respect to <f>. 

In this work, in order to proceed further, we simply treat two special cases for the potential 
function, namely when there is no potential and when there is a non-zero constant value potential, 
V = V Q . 

In free potential case, solutions to equations (7), corresponding to the values of index curvature, 
with the Hamiltonian constraint, Ti, ~ 0, are as follows 



k = 1 : { 



a(t) = Ai cost + A2 sin t and x(t) = B\ cost + B2 sint , 
with constraint : A\ + A\ = B\ + B\ , 



(8) 



k = -1 : 



a(t) = A^e 1 + A 4 e"* and x(*) = 5 3 e * + 5 4 e~* , 
with constraint: A3A4 = B3B4 , 



(9) 



k = 0: { 



o(t) = A 5 t + and x(*) = B ^ + B & 

with constraint : A 2 = B 2 , 



(10) 



where ^Vs and B^s are constants of integration. 

We will compare these solutions with their NC analogues in the next section, where we will also 
discuss the case of non-zero constant potential along with its NC correspondent. 



2.2 Noncommutative Phase Space 

Noncommutativity in classical physics is described by the Moyal product law (shown by the * 
notation in below) between two arbitrary functions of phase space variables, namely £ a = (x 1 ,^) for 
i = !,■■■ ,1 and j = I + 1, • • • , 21, as [7] 



(f*9)(0 = exp 



1 



-a 



ab 



/(CiMc 2 ; 



Ci=C2=C 



(11) 
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such that 

(«-*)=(_/* Sij + aij ), (12) 

where a,b = 1, 2, • • • , 21, 0ij and fyj are assumed to be elements of real and antisymmetric matrices, 
G{j is a symmetric matrix (which can be written as a combination of 0^ and Pij), and dimension of 
the classical phase space is 21. The deformed or modified Poisson brackets are defined as 

{f,g}* = f*g-g*f, (13) 

and also the modified Poisson brackets of variables are 

{xi, Xj} a = 9ij, {xi,pj} a = Sij + Oij and {pi,Pj} a = 0ij . (14) 

The simplest way to study physical theories within the NCG is replacement of the Moyal product 
with the ordinary multiplication. Actually, where variables of classical phase space obey the usual 
Poisson brackets, i.e. {x{, Xj} = = {pi,Pj} and {x{,pj} = 5{j, one can consider the following suitable 
linear, non-canonical, transformation 

x'i = Xi- -0ijp> and p' i =p i + ^fyx 3 , 

with constants 6^ and , for which the usual Poisson brackets of the primed variables yield 

{x'i, x'j} = dij, {x'^p'j} = Sij + Oij and {p^p'j} = foj , (16) 

where a, L j = —0k(iPj)i8 kl /4. Hence, relations (16) give the same results as relations (14). That is, 
one can still employ the commutative variables instead of the NC ones, where it is more convenient 
to work with the Poisson brackets (16) than the modified Poisson brackets (14). The transformation 
(15) can be considered as an extension of the classical mechanics to the NC classical mechanics. 
In geometrical language, the usual Poisson brackets are mapped to the modified Poisson brackets 
through transformation (15) and it should be noted that the two sets of deformed and ordinary 
Poisson brackets represented by relations (14) and (16) must be considered as distinct relations. 
Although, for a compatible extension, the transformation (15) must have an inverse and this imposes 
some conditions on the NC parameters, which we will specify for our model in below. 

Furthermore, let TL = TL(xi,pi) be the Hamiltonian of a system in the commutative case, we shift 
the canonical variables through (15) and assume that the functional form of the Hamiltonian in the 
NC case is still the same as the commutative one, i.e. 

H nc = W&.pJ) =n(xi- X -Q l3 pP, pi + ■ (17) 

This function is defined on the commutative space and therefore, the equations of motion for unprimed 
variables are obviously x l = dTC nc /dpi and p % = —dTt nc /dxi. Evidently, the effects due to the 
noncommutativity arise by terms including the parameters Oij and (3ij. 

Let us now proceed to study the behavior of the model in a phase space with deformed Poisson 
brackets (16) such that the minisuperspace variables do not commute with each other. For our model, 
the two gravitational degrees of freedom are the minisuperspace coordinates (a, x)- The corresponding 
phase space variables are (x 1 , x 2 , p 1 , p 2 ) = (a, X^Va^Vx)-, where p a and p x are the linear minisuperspace 
momenta. We assume the NC parameters 

6 12 = 9 > and /3 12 = 4/3 > , 

hence a 12 = 6(3, where and (5 are constants. Thus, after making the transformations 



a -> a - -p x , p a ^Pa + Wx , X^X + -jPa and p x -> p x - 2(3a, (19) 
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(15) 



(18) 



in (6), the NC Hamiltonian will be 



1 — k9 2 1 
—(Pi ~ p\) + (kO - (3)( X Pa + ap x ) + ((3 2 - k)(a 2 - X 2 ) + 



where V = V{4>) is a modified potential with <f> = y/2{x + 0p a /2)/{a — 8p x /2). Besides, the inverse 
transformation of (19) exists when its determinant is not zero, that is when 9(3 / 1. Now, the 
equations of motion for the NC Hamiltonian are 

a = ~(i - ke 2 ) Pa + (ko - 0) x + ^=o(a - l -9 Px fv' , 

p a = 2(k - (3 2 )a + {f3- k6) Px - 4(a - \e Px fV + y/2( x + \o Pa ){a - ^9p x ) 2 V' , 

X=(k0- 0)a + 1(1 - k0 2 ) Px - 29(a - \dp x fV + ±=0{x + \op a )(a - \9p x ) 2 V' , 

p x = ((3- k9) Pa + 2{(3 2 - k) X - V2(a - \dp x fV' , (21) 

where V' = dV/d(f>. Comparing equations (21) with (7) shows that, in general, equations of motion 
are coupled in the NC case and for 9 = = (3, equations (21) reduce to (7) as expected. 
In free potential case, equations (21) read 

a = -^(l-k9 2 )p a + (k8-p)x, 
p a = 2(k - (3 2 )a + {(3- k9)p x , 
x = (k8-(3)a + ^(l-k9 2 )p x , 

Px = {f3- k9) Pa + 2{f3 2 - k) X ■ (22) 
By eliminating momenta variables, one gets 

a = -m 2 a + 2(k8 - (3)x and X = ~m 2 x + 2(k0 - (3)a , (23) 
with the Hamiltonian constraint 

(P -1) (a 2 - x 2 ) = constant if 1 - k9 2 = } 

or 

m 2 (a 2 -x 2 )- (a 2 -X 2 ) = if 1 - k9 2 ^ , (24) 

where m 2 = k(l — (39) 2 . It is noticeable that when k = 0, the NC parameter 9 is removed from 
equations (23), that is, in the spatially flat FRW universe, the motion equations are affected only by 
the NC parameter (3. 

If one takes A = (1 — k9 2 ){(3 2 — k), real solutions of equations (23) can be written as 



A > : < 



a(t) = e^-^*(^i sinh y/At + B l cosh y/At) + e'^-^id sinh y/At + D 1 cosh y/At) 



x (t) = e (ke-p)t ( Al sinh ^ + Bl cosh _ e -(M-P)t(d sinh y/At + Di cosh y/At) 

(25) 

and 



A < 



i(t) = e^ e -^\A 2 sin x/^At + B 2 cos y/^At) + e^* 9- ^*^ sin y/^At + D 2 cos y/^At) 



x (t) = e ( k9 -^ t (A 2 sin y/^At + B 2 cos v /Z Ai) - e^ ke -^\C 2 sin + D 2 cos y/^At), 

(26) 

where Aj's, B^s, Q's and Dj's are constants of integration. 



1 Note that, the condition 1 — k6 2 — is possible only when k — 1 and hence 6 = 1. 
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The quantity A is always positive when k = 0,-1, hence, in general, solutions (25) are in the 
form of 

2 

]T (beW + ce-**) , (27) 

3=1 

with b and c as new constants, and 

6i = k9 - (3 + and S 2 = k9 - (3 - VA . 

Comparing solutions (27) with its commutative analogues (9) and (10) shows that for k = — 1, 
solutions are still hyperbolic, but with an extra coefficient in the exponent that depends on the NC 
parameters and gives enough room for better adjustments. In the case k = 0, one has 8± = and 
82 = —2/3, and hence, the time dependence of solutions have been modified from linear in commutative 
case to hyperbolic in the NC case, where the former solution is not suitable for an accelerating universe, 
though the latter one is capable to be adjusted with the observed accelerated expansion. Also, in the 
late time, solution (27) describes a de Sitter universe for which its cosmological constant is written 
in terms of the NC parameter as A/3 = A(3 2 . 

For k = 1 geometry, one has A = (1— 9 2 ){0 2 — 1), which can be positive or negative or zero, depends 
on the different choices of the NC parameters. That is, when (0 > 1, (3 < 1) or (9 < 1, (3 > 1), A 
is positive. Hence, we have solutions (25), however, in this case the NC parameters have upper and 
lower bounds. When A < 0, existence of oscillating solutions is provided when both NC parameters 
simultaneously have a lower bound, namely 9 and (3 > 1, or an upper bound, namely 9 and (3 < 1. 
The commutative solutions (8) are oscillating with the period of 2tt and constant amplitudes, whereas 
the NC solutions (26) have the period of 2tt/\/—A and varying amplitudes with factor e^ - ^*. Note 
that, the choice 9 = (3, leads to a constant amplitude as the commutative case, but with a period 
of 27r/ 1 1 — 9 2 \. The case A = is possible when 9 = 1 or when (3=1, where solutions are again 
hyperbolic. The case k = 1 and 9 = 1, resembles the condition 1 — k9 2 = with arbitrary (3, however, 
with the value of (3 = 1 one gets trivial constant solutions. 

Now, we treat a non-zero positive constant potential V = V . Since, in general, the motion 
equations are not sufficiently simple to be solved, we restrict ourself to the case k = and (3 = 0. 
Thus, after a little algebra, equations (21) give 

d 2 = V D {a - v Q 9f + d and X {t) = Xo + v Q t - 9a , (28) 

where Xo, d and v Q = p x /2 are constants and the Hamiltonian constraint gives d = v 2 . 
By setting = 0, one obtains the commutative equations as 

a 2 = V Q a 4 + v 2 Q and x (t) = X o + v Q t . (29) 

Solutions of a(t) are in terms of the Jacobi elliptic functions, i.e. 

a(t) = ± D sn (iy/VoDt, i) , (30) 

with D = [v / (i\/Vo )] x / 2 and the initial condition a(0) = 0. These two solutions are periodic functions 
with respect to the time, and for Vo < 0, they qualitatively behave as sine functions with amplitude 
\D\. For Vo > 0, solution (30) qualitatively behaves as a tangent function and its diagram is plotted 
in Fig. 1 (left) for numerical values v Q = 1, Xo = 2, = 2 and V Q = 1. Obviously, physical solution 
corresponds to the positive part of this diagram. 

Equations (28) and (29) show that the effect due to the noncommutativity is a shift in the scale 
factor with a magnitude \ v Q 9\ that changes the zero point of the scale factor. This effect can remove a 
negative scale factor when the constant potential is negative. That is, in such a case, the corresponding 
NC solution is vq9 ± D sn {i\f\T Dt, %) , where by choosing v$9 > \D\ with vq > 0, then the resulted 
scale factor will be positive. Clearly, plot of a(t) in the NC case is the same as the commutative one 
with a vertical transfer of magnitude. 
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The linear behavior of the scalar field in the commutative case is altered by an additional term 
— 6a. Two plots of the scalar field in the NC case are sketched in Fig. 1 (right) for numerical values 
v = 1, Xo = 2 and 9 = 2, one with V$ = — 1 (solid line) and one with Vq = +1 (dashed line). Note 
that, in a positive potential, the scalar field goes to ±oo for very late time (when t goes to infinity), 
whereas, in a negative potential, it goes to +oo monotonically. 

In terms of the cosmic time, dr = adt, the first equation in (29) may have a solution, for Vq > 0, 

as 

a(r) oc [ sinh {2y/V Tj ] 1/2 . (31) 
The behavior of such a scale factor in the late cosmic time, i.e. r S> 0, is 



a(t) oc 



U -t 



(32) 



where t\ is a constant such that if r — > oo, then £ — > tj. Hence, in this limit, the scalar field in the 
NC case grows proportion to 6 as 6(ti — t)~ 2 , and deviation from the linear time dependence becomes 
very large, as one can recognize it in Fig. 1 (left). 



3 The Quantum Model 

We proceed to quantize the cosmological model given by the action (1) in the case of free potential, 
such that the canonical quantization of the phase space leads to the Wheeler-DeWitt (WD) equation, 
"H^ = 0, where H is the Hamiltonian operator and ^ is a wave function of universe. For arguments 
about the quantization based on WD equation in the FRW universe including matter fields, see, 
e.g., Refs. [15]. We employ the usual canonical transition from classical to quantum mechanics 
via the generalized Dirac quantization from the Poisson brackets to the quantum commutators, i.e. 
{} — > Then, as the classical approach, we investigate the commutative and NC frames in the 

following subsections. 



3.1 Commutative Frame 



As usual, the operator form of Hamiltonian (6) can be acquired by the replacements p a — ► — id a 
and p x — ► — id x . Assuming a particular factor ordering, the corresponding WD equation, for V = 0, 
is 

^ <9 2 ... 2 2 ; 



*(a, X ) =0. 



Considering the following change of variables 

a = p cosh ip and X = P smn V > 

equation (33) reads 

(&_ + ld_ _ 1_^ _ Akp2 \ 

I dp 2 p dp p 2 d(p 2 J ' 

We assume solutions of equation (35) as a product ansatz 

V(p,<p)=iP(p)e 2 ^, 
where a is a real constant. Substitution of (36) in equation (35) leads to 

p V/r / 



(33) 

(34) 
(35) 

(36) 
(37) 
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where the prime denotes derivative with respect to p. The well-defined eigenfunctions of equation 
(37), considering boundary conditions, can be written as 



Jia{-kp 2 ) 

cos(2aln p) 
{ K ia {kp 2 ) 



for k = — 1 
for k = 
for k = 1 , 



(38) 



where the functions J u and K v are the first and modified second kind of the Bessel functions, respec- 
tively. Hence, the wave packet corresponding to (38) is 



/+oo 
C a Mp)e 2iaip da, 
-oo 



(39) 



where C a can be taken [9, 10] to be a shifted Gaussian weight function with constants b and c, i.e. 

e -b(a~c)\ 

Fig. 2 shows plots of the probability |^| 2 for k = and k = 1 with values 6 = 1 = c, in the range 
of < p < 10 and — 5 < ip < 5. The corresponding plot for k = — 1 is qualitatively similar to the case 
k = 1. 



3.2 Noncommutative Frame 

We assume that, in a general NC quantum phase space, the coordinate operators of the FRW 
minisuperspace and their generalized momenta obey the star deformed Heisenberg algebra, like the 
ones in noncommutative quantum mechanics, as [16] 

Xj] a = i9ij, [xi,pj] a = i(5ij + Gij) and \pi,Pj] a = ifaj ■ (40) 

The notations and definitions are the same as in the NC classical model. Now, the corresponding 
noncommutative WD equation can be written by replacing operator product with the Moyal product, 
namely 

Hnc(xi,p i )*V(x,p) = 0. (41) 

It is well-known in noncommutative quantum mechanics [16] that the original phase space and its 
symplectic structure can be modified to reformulate relations in the commutative algebra when the 
new variables, as in (15), are introduced. Hence, the original equation (41) for these new variables 
reads [17] 

n QC (xi - \ei$,Pi + \Pij&) v(x, P ) = o. (42) 

Therefore, the noncommutative WD equation corresponding to relation (20), with 6(3 ^ I and for 
V = 0, can be written as 

"(1 - k6 2 )(d 2 a - d 2 ) - Ai{k6 - [3)( X d a + ad x ) + 4(P 2 - k){a 2 - X 2 ) 

that, with the change of variables (34), reads 

<' - k ° 2 i$ + M - 7w) - ^ - < - 4 < fc - ^ 

If 1 — k6 2 ^ 0, by the ansatz (36), equation (44) reduces to 

V>" + — + 4 ( ^ - ep 2 + 2 7 a J V = , (45) 



*(a, X )=0, (43) 



* = . (44) 
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where e = (k — (3 2 )/(l — k9 2 ) and 7 = (k6 — [3)/{l — kd 2 ). Before discussing solutions of (45), we 
enumerate a few features of this equation. Firstly, it reduces to equation (37) when one sets e = k 
and 7 = 0, where a trivial case is 9 = = (3, as expected. Secondly, it is interesting that, in 
the two following cases where the noncommutativity is still present, equation (45) again reduces to 
equation (37). These two cases 2 are k = 1 with 9 = (3 ^ 1, and k = with (3 = and arbitrary 
9. In another words, in these cases, general form of the NC wave functions are the same as the 
commutative solutions (38). Finally, in the case k = 0, solutions to equation (45) does not depend 
on the 9 parameter, which is also a common feature in the classical model. 

A particular solution of equation (45) can be written in terms of the Whittaker functions, ;U 
and M^, as 

^ (p) = P' 1 [i^ ,u (2vV) + B a W„ ,„ (2vV) ] , (46) 

where A a and B a are superposition constants, p = 07/ ^fl and v = ia. The both Whittaker functions, 
even in classically forbidden regions, are convergent when e is negative^ that can obviously occur for 
k = and k = — 1. Also, e can be negative in the case k = 1 for special values of the NC parameters, 
we will discuss these situations in below. 

In the case k = for which e = — /3 2 , one can write solution (46) with one of the Whittaker 
functions, e.g., M^^ and gets its corresponding wave packet as 

(p, <p) = p- 1 r e~ b ^ 2 M ia , ia (2z0p 2 ) e 2 ^da . (47) 

Fig. 3 shows the probability |^| 2 corresponding to (47) for (3 = 1/10 and the Gaussian constants 
b = 1 = c. The comparison of this figure with the left diagram in Fig. 2 identifies that, although 
the most probable state of a commutative universe is around ip = 0, the most probable state related 
to a NC universe is shifted for ip less than zero. Also, another difference in those two figures is the 
displacement of peaks to greater values of p in a NC universe with respect to a commutative one. 
By consecutive diagrams, it can be shown that this difference is also more manifested when (3 gets 
smaller values. 

As mentioned before, the case e < can also be attained when k = 1, for which e = (1— j3 2 )/{l— 6 2 ). 
This quantity is negative when one of the inequalities (9 < 1 and (3 > 1) or (0 > 1 and (3 < 1) are held. 
The first index of the Whittaker functions, for both of these inequalities, is also a pure imaginary 
number. When Q{3 = 1, wave packets resemble those of the case k = 0, but this choice violates the 
existence of the inverse transformation for (19). On the other hand, plots of wave packets for 9/3^1, 
in general, show that they are not suitable for description of a real universe. However, choosing special 
values, such as considering the first inequality with 9 = and (3 very large values that [3 2 ± 1 rs (3 2 
or considering the second inequality with [3 = and 9 very large values that 9 2 ± 1 « 9 2 , one can 
get p ~ ia and hence again, construct well-behaved wave packets in a closed universe, which their 
corresponding plots are qualitatively similar to Fig. 3. For instance, Fig. 3 does also illustrate a 
similar plot of probability for (3 « and 9 « 10 in this case, where the visible contrast between such a 
figure and its commutative analogue, Fig. 2 (right), implies that the probability of expansion in a NC 
universe is more than a commutative counterpart. Also, the most probable state of a commutative 
universe is around <p = and p = 0, but the most probable state related to a NC universe is shifted 
for <p less than zero and p greater than zero. Besides, there is displacement of peaks to greater values 
of p in a NC universe with respect to a commutative one. This effect is more manifested when 9 gets 
greater values. Hence, there are different possible universes (states) from which, our present universe 
could be evolved and tunneled in the past from one state to another^ 

The above discussion, in the latter paragraph, is completely fulfilled when k = — 1, except that 
there is no bound on the NC parameters, for e is always negative in this case. 



2 The case k = — 1 with 9 — —(3 is not valid, for our assumption (18). 

3 For properties of the Whittaker functions and their indices, see, e.g., Refs. [18]. 

4 Such a kind of arguments can also be found in, e.g., Ref. [10]. 
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For when e is positive, the Whittaker term of M„ v , in classically forbidden region, is divergent. 
Besides, the both Whittaker terms, when e > 0, are proportional to exp(— y/ep 2 ), and hence, are 
quickly damped as p grows. 

We should mention that e also vanishes for the case k = 1 with (3 = 1, which leads to 7 = 
— 1/(1 + 6). In this case, solution of equation (45) can be considered as 

ip a (p) (x J ia (2ijx/ap) . (48) 

This solution, as can easily be checked, is not well defined for constructing wave packets and hence, 
we do not consider it. 

The condition 1 — k9 2 = is equivalent to k = 1 and 9 = 1, for which equation (44) reduces to 



i(l-/3)^ + (l-/3V 



* = . (49) 



Its solution is 

* = R{p)e i{l+ ^P 2 , (50) 

where R(p) is any differentiable function. The probability of this wave function is independent of the 
(3 parameter. 



4 Conclusions 

In this work, we have carried out an investigation into the role of NCG in cosmological scenario 
by introducing a NC deformation in the algebra of phase space variables almost along the same lines 
proposed in Ref. [12]. The space is generated by a conformal scalar field when is non-minimally 
coupled to geometry whose background is the FRW metric, in all three cases k = —1,0,1. The 
noncommutativity has been introduced in the both, space and momentum, sectors with more attention 
on the role of the NC parameter (3. The investigation has been carried out by means of a comparative 
analysis of the evolution of universe and the most probable universe wave functions in classically 
commutative and NC frames and quantum counterparts. 

In absence of the scalar field potential function, the classical exact solutions have been obtained 
in both commutative and NC frames. Contrary to the commutative case, the NC solutions can be 
regulated with the aid of both NC parameters. Especially, in spatially flat universes, the major player 
in tuning solutions is the (3 parameter, where it can be employed to adjust time dependent solutions 
with the observational data, e.g., the accelerating expansion rate. Interestingly, we have found that 
the existence of particular solutions imposes bounds on the values of NC parameters. 

In the presence of a non-zero constant potential and absence of f3, the time dependence of the 
scalar field changes contrary to the scale factor which is just shifted by a constant value proportional 
to the 9 parameter. This shift can remove negative scale factors when negative constant potentials are 
used. The scalar field deviation from linear time dependence of the commutative case, for a positive 
constant potential, becomes more apparent in the late cosmic time. 

One expects that when noncommutativity effects are turned on in the quantum scenario they 
should introduce significant modifications. To see this, the corresponding quantum cosmology has 
also been considered and the exact solutions, through the WD equation in the commutative and NC 
frames, are obtained in order to investigate differences in the solutions. As in the case of classical 
model, the (3 parameter is the only one that is responsible for NC effects in spatially flat universes. 
Comparing the numerical diagrams of the probability in commutative and NC cases, especially in 
the case of positive and negative curvature indexes, shows that expansion in a NC universe is more 
probable, depending on the values of the NC parameters, than a commutative counterpart. In the 
flat case, existence of universes with more possible states invokes smaller values of (3 in NC frames. 
Also, existence of more probable states and transferring peaks of probability to greater values of p, 
defined in (34), in the NC frames with respect to the commutative ones are intensified when (3 gets 
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smaller values. Also, the latter property can qualitatively be apparent in the closed and open cases 
when, for instance, (3 = with sufficiently large values, or 6 = with (3 sufficiently large values. 
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Figure 2: The probability \^/\ 2 in the commutative case for k = (left) and k — 1 (right), both 
with the Gaussian constants b = 1 — c. 




